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Abstract 

We present a few combinatorial identities which were encountered in our work 
on the spectral theory of quantum graphs. They establish a new connection 
between the theory of random matrix ensembles and combinatorics. 

I. INTRODUCTION 

In the present paper we show that some questions arising in the study of spectral correlations 
for quantum graphs can be cast as combinatorial problems. Solving these problems for a 
particular system, we discovered the following novel combinatorial identities: 

(i) Let n, q be arbitrary integers with 1 < q < n and 

^-)-^^;,H::i)(r-y"- , - 1 ^"- , - lv 

Then 



min(g,n— q) 

S{n,q)= F uy (n,q) = l. (1) 



(ii) Let s,t be arbitrary positive integers and 

t) = m E\-±y- u (*) (* : j) = (-i) s+< ( s + 1 ~ x ) 1,2 p s+t -i,s(t) , (2) 

where -P/v^x) are the Kravtchouk polynomials to be defined in Eq. (f46[). Further, let 
x,y be complex with \y\ < l/y/2. Then we have the generating functions 

Gl( ,) = ±* M *« = ^ - ^) , (3, 



G 2 (x)= Y Af(s, t) Af(t, s) x s+t =- + + - - (4) 
W ^ 1 j 2 (2x + 1)^4x2 + 1 2 i; 



and 



(iii) Let m be any positive integer. Then 



q=l 

and 



w E |«' =2 »* + (-.rh-2 (6 ) 



(2m + lff;(^Y^J =2 2m+2 -2(-l)™^ -2. (7) 



(iv) Let < g < n, and define 



Afri ) = 1 / 1 for g = 0, n , . 

\/2" I (-l) g (n/q) Af(n - q,q) for < g < n . ^ 



Then for any positive integers < k < ^ and an arbitrary integer no, 

oo n 

lime y e~ ne yA(n + v,q + n)A(n,q) = A(v,k) . (9) 

These identities establish a new connection between combinatorics and the theory of random 
ensembles of 2 x 2 matrices. The physical background and motivations are described in a few 
recent publications , to which the interested reader is referred. An immediate application 
of ([!]) is also given in @. Here, we shall provide the minimum background necessary for 
the understanding of the combinatorial aspects of the problem, and for a self-contained 
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exposition of our results. It will also enable us to propose a conjecture which generalizes the 
combinatorial approach to random matrix theory for matrices of large dimensions. 

We start with a few definitions: Graphs consist of V vertices connected by B bonds (or 
edges). The valency Vi of a vertex i is the number of bonds meeting at that vertex. We shall 
assume that two vertices can be connected by a single bond at most. We denote the bonds 
connecting the vertices i and j by b = The notation will be used whenever we do 
not need to specify the direction on the bond. Hence [i,j] = [j,i]- Directed bonds will be 
denoted by d = and we shall always use the convention that the bond is directed from 

the first to the second index. If d = (i,j) we use the notation d = (j, i). Let gw be the set 
of directed bonds (i, j) which emanate from the vertex i, and gW the set of directed bonds 
(j', i) which converge at i. The vertices i and j are connected if g^' H g^) 7^ 0. The bond d! 
is connected to the bond d if there is some vertex i with d e gW and d! E g^ 1 '. d and d are 
always connected. 

The Schrodinger operator on the graph is defined after the natural metric is assigned 
to the bonds, and the solutions of the one- dimensional Schrodinger equation on each bond 
(id x — A) 2 tp(x) = k 2 ifj(x) is given as a linear combination of counter-propagating waves. 
A stands here for a magnetic flux. The (complex) amplitudes of the counter propagating 
waves are denoted by a^, where the subscript d stands for the directed bond along which 
the wave propagates, d = 1,..,2B. Appropriate boundary conditions at the vertices are 
imposed, and the spectrum of the Schrodinger operator on the graph is determined as the 
(infinite, discrete) set of energies k 2 , for which there exists a non-trivial set of ad which is 
consistent with the boundary conditions. The condition of consistency can be expressed by 
the requirement that 

det(J - S(k)) = (10) 

where the bond- scattering matrix S(k) is a unitary operator in the Hilbert space of 2B 
dimensional vectors of coefficients a^. The unitarity of S(k) ensures that the spectrum of 
the Schrodinger operator is real. The matrix S(k), which is the object of our study is defined 
as 

^(«),(J,m)(*)=^ iJ ' ( *Vg:(fc). (11) 

The matrix elements Sd,d'(k) vanish if the bonds are not connected. As a consequence, 
the unitarity of S implies also the unitarity of the fj-dimensional vertex-scattering matrices 
cr}^ and vice versa. The phases (fiuj), are given in terms of the bond length Luji, and the 
magnetic flux A {iJ) = -A {jti) QTJ , 

<%i)( fc ) = ( k + A (i,j)) L liJ] ■ ( 12 ) 

The two phases pertaining to the same bond <pd and <p d are equal when Ad = A d = 0. In 
this case S is symmetric and the Schrodinger operator on the graph is invariant under time 
reversal. Time-reversal symmetry is violated when some magnetic fluxes do not vanish. 

S can also be interpreted as a quantum time evolution operator describing the scattering 
of waves with wave number k between connected bonds. The wave gains the phase (j>(i,j){k) 
during the propagation along the bond while the a^m describe the scattering at the 
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vertices. In this picture, the unitarity of S guarantees the conservation of probability during 
the time evolution. 

We will avoid unnecessary technical difficulties and consider the matrices af^ to be 

k independent constants. One may find explicit expressions for crj^i by requiring besides 
unitarity that the wave function is continuous at the vertices. The resulting expression is |I| 

2 

crj*]/ = 8j t f (Neumann b. c.) . (13) 

Note that back-scattering (j = j') is singled out both in sign and in magnitude. In all 
nontrivial cases (v j > 2) the back-scattering amplitude is negative, and the back-scattering 
probability |cr^]| 2 approaches 1 as the valency Vi increases. 

Finally, a "classical analogue" of the quantum dynamics can be defined as a random 
walk on the directed bonds, in which the transition probability between bonds 
connected at vertex j is Icp^l 2 . The resulting classical evolution operator with matrix 
elements 

U d4 , = \S d>d ,\ 2 (14) 
is probability conserving, since unitarity implies J2i Wul 2 = 1- 



II. SPECTRAL TWO-POINT CORRELATIONS AND PERIODIC-ORBIT SUMS 

The spectrum of S consists of 2B eigenvalues e tBl which are confined to the unit circle. 
Their distribution is given in terms of the spectral density 

2B nn -i oo 

d{6) = - ft) = t- + t- E s ^~ Wn + c - c - > ( 15 ) 
1=1 2n 2n „=1 

where 8^ denotes the 2ir periodic delta function. The first term on the r.h.s. is the average 
density d = §f ■ The coefficients of the oscillatory part 

s n = trS n 

will play an important role in the following. s n = trS n is a sum over products of n matrix 
elements of S, and because of (O) the bond indices of each summand describe a connected 
n-cycle (= n-periodic orbit) on the graph 

s n ='£A p j*>. (16) 

pev n 

In C p~6|) V n denotes the set of all n-periodic orbits (PO's) on the graph. Note that for the 
convenience of presentation we will consider cycles differing only by a cyclic permutation as 
different PO's. The phases $ p = X^=d 4>dj can be interpreted as the action along the PO p. 
The amplitudes A p are given by 

n-l 

A=T[s dj+udj , (i7) 

3=0 
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where j is understood mod n. Sometimes it is useful to split the amplitude in its absolute 
value and a phase e 1Mp7r . For example, in the case of Neumann b. c. (13) \i v is an integer 
counting the number of back-scatterings along p. 

In complete analogy to (p~6|) we can represent also the traces of powers of the classical 
evolution operator 



u n = trU r 



(18) 



as sums of periodic orbits of the graph. 

The two-point correlations in the spectrum of S flT5|) can be expressed in terms of the 
average excess probability density R2(r; (3) of finding two phases at a distance r, where r is 
measured in units of the mean spacing 



R2(r; 0) = (^ f_l ^ d (6) d(9- [tt/B] r] 



B 



2 - / 7T 

— cos I —nr 

^ n n=l 



B 



1 

2B 



(19) 



The bond scattering matrix depends parametrically on the phases <pd (0). We shall define 
two statistical ensembles for S in the following way. The ensemble for which time-reversal 
symmetry is broken consists of S matrices for which the cbd are all different, and we consider 
them as independent variables distributed uniformly on the 2B torus. Invariance under time 
reversal implies (f>d = <Pd anc ^ ^ ne corresponding ensemble is defined in terms of B independent 
and uniformly distributed phases. We shall distinguish between these ensembles by the value 
of the parameter (3 = {number of independent phases}/!?. Expectation values with respect 
to these measures are denoted in (|T9| ) by triangular brackets, 



0B 

(•••)^n 



i 

2^ 



(20) 



The Fourier transform of R2{r; (3) is the form factor 



K{n/2B;(3) 



2B 



(21) 



on which our interest will be focussed. If the eigenvalues of the S were statistically indepen- 
dent and uniformly distributed on the unit circle, K(n/2B) = 1 for all n. Any deviation of 
the form factor from unity implies spectral correlations. Using (|16D the form factor (|2lD is 
expressed CIS db double sum over PO's 



K(n/2B;f3) = — 
1 

" 2B 



E A 



peVn 



P,p'ev„ 



(22) 



P 



In order to perform the average over all the phases <pd in (|22|) we write 
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$p = E4 P) ^, (23) 
d 

where rffi counts the number of traversals of each directed bond such that Y,d n ^ = n - 
According to ( p0|) we have 

e i(n0 d +n'^)^ =i = 6 nfi 6 n , fi , (24) 

; iWi ^^ =2 = <W«,o- (25) 

Thus, the double sum in (53) can be restricted to families of orbits. For [3 = 2, let £ 

(C) 

be the family of isometric PO's which have the same integers n d . That is, the family 
consists of all the PO's which traverse the same directed bonds the same number of times, 
but not necessarily in the same order. In the case (3 = 1, £ consists of all PO's sharing 
nl = njj + . That is, the family contains all PO's which traverse the same set of 
undirected bonds the same number of times, irrespective of direction or order. We find 

K(n/2B;(3) = ±- £ |£AI 2 - (26) 
zn c&r{ff) n pec 

J~{(3) n denotes the set of all vectors £ = [rid] for (3 = 2 (£ = [n b ] for (3 = 1) of (3 B non- 
negative integers summing to n, for which at least one PO exists. For Neumann b. c. (|13|). 
e. g., (|26D amounts to counting the PO's in a given set £ = [rid] taking into account the 
number of back-scatterings along the orbit. The problem of spectral statistics is now reduced 
to a counting (combinatorial) problem which is, however, very complicated in general. Even 
the determination of the number of families L for a given n is difficult. For (3 = 2 an obvious 
necessary condition for the existence of a PO with a given set of bond traversals £ = [rid] is 
that at any vertex the number of incoming and outgoing bonds is the same, i. e. 

Yl n d = J2 U d (i = i, . . . V) . 

For (3 = 1, the analogous condition reads 

n d~ n d mod 2 = (i = 1, . . . V) , 

i. e. the total number of traversals of adjacent bonds should be even at each vertex. However, 
it is not so easy to formulate a sufficient condition for the existence of a PO given a set of 
numbers rid- In particular one must take care to exclude cases, in which the set of traversed 
bonds is a union of two or more disconnected groups ("composite orbits"). 

Extensive numerical work II] revealed that for fully connected graphs (vj = V — 1), and 



for V 1, the form-factor (21) is well reproduced by the predictions of random matrix 
theory for the Circular Orthogonal Ensemble (COE) {(3 = 1) or the Circular Unitary 
Ensemble (CUE) {(3 = 2). This leads us to expect that ( ^6|) approaches the corresponding 
random matrix prediction in the limit V — > oo. This conjecture is proposed as a challenge 
to asymptotic combinatorial theory. 
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III. THE RING GRAPH 



In the following we will evaluate explicitly the quantities introduced in the previous 
section for one of the simplest quantum graphs. It consists of a single vertex on a loop 
(see fig. |l|). There are two directed bonds d — 1 and d = 2 with <pi ^ <p 2 , i- e. time-reversal 
symmetry is broken. Since this graph would be trivial for Neumann b.c. the vertex-scattering 
matrix at the only vertex is chosen as 

, N / cos 77 ishi7 
a(r]) = . . 

V 1 Sill 7] COS 1] 

with < rj < 7r/2. The corresponding bond-scattering matrix is 

S( ) = ( ^ ° ^ ( cosr l isin7]\ 
I e^ 2 I I i sin 77 cos f] J 

We shall compute the form factor for two ensembles. The first is defined by a fixed value of 
77 = 7r/4, and the only average is over the phases (ftd according to (f24|). The second ensemble 
includes an additional averaging over the parameter rj. We will show that the measure for 
the integration over 77 can be chosen such that the model yields exactly the CUE form factor 
for 2x2 random matrices ||. 




A. Periodic Orbit Representation of u n 

We will first illustrate our method of deriving combinatorial results from the ring graph 
in a case where a known identity is obtained. Consider the classical evolution operator U of 
the ring graph. According to flH| ) we have 

/ cos 2 77 sin 2 77 \ 
V sin 77 cos T] J 

The spectrum of U consists of {1, cos 27/}, such that 

Un (r]) = 1 + cos"27/. (30) 

We will now show how this result can be obtained from a sum over the periodic orbits of the 
system, grouped into families of orbits as in (P6f). In the classical calculation it is actually 
not necessary to take the families into account, but we would like to stress the analogy to 
the quantum case considered below. The periodic orbit expansion of the classical return 
probability can easily be obtained from Q2D| ) by expanding all matrix products in flTBp. We 
find 

«n= E •■■ E fiiWifa). ( 31 ) 

ii=l,2 i„=l,2 j=0 

where j is again taken mod n. In the following the binary sequence [ij] (ij G {1,2}; j = 
0, . . . , n — 1) is referred to as the code of the orbit. We will now sort the terms in the multiple 
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sum above into families of isometric orbits. In the present case a family is completely 
specified by the integer q = q\ which counts the traversals of the loop 1, i.e., the number of 
letters 1 in the code word. Each of these q letters is followed by an uninterrupted sequence 
of tj > letters 2 with the restriction that the total number of letters 2 is given by 

Y i tj = n-q. (32) 
i=i 

We conclude that each code word in a family < q < n which starts with i\ = 1 corresponds 
to an ordered partition of the number n — q into q non-negative integers, while the words 
starting with i\ = 2 can be viewed as partition of q into n — q summands. 

To make this step very clear, consider the following example: All code words of length 
n = 5 in the family q = 2 are 11222, 12122, 12212, 12221 and 22211, 22121, 21221, 22112, 
21212, 21122. The first four words correspond to the partitions + 3 = 1 + 2 = 2 + 1= 3 + 
of n — q = 3 into q = 2 terms, while the remaining 5 words correspond to 2 = + + 2 = 
+ 1 + 1 = 1 + + 1 = + 2 + = 1 + 1 + = 2 + + 0. 

In the multiple products in (pTF) a backward scattering along the orbit is expressed by 



two different consecutive symbols ij ^ in the code and leads to a factor sin 2 ?], while 
a forward scattering contributes a factor cos 2 n . Since the sum is over periodic orbits, the 
number of back scatterings is always even and we denote it with 2u. It is then easy to see 
that v corresponds to the number of positive terms in the partitions introduced above, since 
each such term corresponds to an uninterrupted sequence of symbols 2 enclosed between two 
symbols 1 or vice versa and thus contributes two back scatterings. For the codes starting 
with a symbol 1 there are (j^j ways to choose the v positive terms in the sum of q terms, 

and there are ways to decompose n — q into v positive summands. After similar 

reasoning for the codes starting with the symbol 2 we find for the periodic orbit expansion 
of the classical return probability 



ra-1 



u n (r ) ) = 2cos 2n r ] + J2H 




q=l v 

n I q — 1\ In — q — V 



sm 4u r] cos 2n ~ 4 ^ 



n— 1 



2 cos 2 ™ n + E E i[ I _ i ( „ ! / ^ cos 2 — r, (33) 



<j=l v 

The summation limits for the variable v are implicit since all terms outside vanish due to 
the properties of the binomial coefficients. From the equivalence between (^) and ( |3"3"D the 
combinatorial identity 

could be deduced which indeed reduces (RBI) to a form 



u n (ri)=2Y:(^ sin 4 ^ cos" 



7/ 

cos 2 r/ + sin 2 77) n + (cos 2 ?7 — sin 2 ?])™ , (35) 



S 



which is obviously equivalent to (pOj). 

( |34|) can also be derived by some straight forward variable substitutions from the identity 



n—m 



, , y/y y m y y/ + m + l 



which is found in the literature 



B. Quantum Mechanics: Spacing Distribution and Form Factor 

In the following two subsections we derive novel combinatorial identities by applying the 
reasoning which led to (|34]) to the quantum evolution operator (^) of the ring graph. We 
can write the eigenvalues of S(r/) as e 1 ^ 1+ * 2 - ) / 2 e ±lA / 2 with 



A = 2 arcos 



COS Tj COS 



(37) 



denoting the difference between the eigenphases. For the two-point correlator we find 



R 2 {r,rj) 



+ 7T 



de d e 



7r r 



die 



nr 



8 2 {r) - 1 / 5 27T (tt r + A) + S 27r r — A) 



7T 



<t>l,2 



5 2 (r) — 1 sin |7rr/2| @(cos 2 r] — cos 2 (7rr/2)), 



7T 



271 



\J cos 2 ri — cos 2 (7rr/2) 



(38) 



Here, 52 (r) is the 2-periodic 5 function. In particular for equal transmission and reflection 
probability (rj = n/A) we have 



R 2 (r,n/A) 



52 (r) 



71 



1 1 



cos (nr) 



cos 7rr 



(39) 



and, by a Fourier transformation, we can compute the form factor 



K(n,n/A) = n dr cos (mxr) R 2 (r,7i/A) 



1 + 
1 + 



-1 



\m+n 



2m 



2 2m+l ^ m 



3 A 

2 n '° 



2\/7vn 



(n > r 



(40) 
(41) 



where m = [n/2] and [■] stands for the integer part. 

Next we consider the ensemble for which transmission and reflection probabilities are 
uniformly distributed between and 1. For the parameter 7] this corresponds to the measure 
dp,(rj) = 2\ cos 7/ shield//. The main reason for this choice is that upon integrating (|38|) one 

gets 
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^av) (r) = Mr)-l + sin (7TT-/2) 

7T 7T 

which coincides with the CUE result for 2x2 matrices. A Fourier transformation results in 
The form factors (|40|) , fl4lf) and ( |43|) are displayed in Fig. [I] below. 



C. Periodic Orbit Expansion of the Form Factor 



An explicit formulation of ([26|) for the ring graph is found by labelling and grouping 
orbits as explained in the derivation of We obtain 



ft 

K 2 (n;r]) = cos 2n r? + — J2 

1 q=l 



2 n-1 



E 



:-i) 1 



V 



'q-r 



n — q — 1 
. v- 1 . 



2f n—lv 

sin cos 77 



(44) 



where g denotes the number of traversals of the ring in positive direction and 2v is the 
number of backward scatterings along the orbit. The inner sum over v can be written in 
terms of Kravtchouk polynomials as 



2 n— 1 

K 2 (n;r]) = cos 2n r/ + - ^ 

1 q=l 



'n-V 
n-q. 



cos 2q r]sm 2{n - q) ri 



n 



-P, 



(cos 2 »),sin rf) 



n—l,n—q 



(45) 



and the Kravtchouk polynomials are defined as in 



by 



Pn± \ x ) 



. k 1 



[uv\ 



1-1/2 k 

u 

u=0 



\k—u 



X 



'N 
. k - 



• x 



u k ~ v v u 



< k < N 

U + V = 1 



(46) 



These functions form a complete system of orthogonal polynomials of integers x with < 
x < N. They have quite diverse applications ranging from the theory of covering codes 
to the statistical mechanics of polymers J7|. The same functions appear also as a building 
block in our periodic-orbit theory of Anderson localization on graphs 0. Unfortunately, 
we were not able to reduce the above expression any further by using the known sum rules 
and asymptotic representations for Kravtchouk polynomials. The main obstacle stems from 
the fact that in our case the three numbers N, k, x in the definition ( |46D are constrained by 
N — k + x — 1. 

We will now consider the special case rj = tt/A for which we obtained in the previous 
subsection the solution (f4"0"p. The result can be expressed in terms of Kravtchouk polynomials 
with u = v = 1/2 which is also the most important case for the applications mentioned above. 
We adopt the common practice to omit the superscript (u, v) in this special case and find 



1 1 

K 2 (n; 7i74) = 1 - 



n-l 

E 

9=1 



fn — 


3 


\n — 





-p 



1 2 



Q 



(?) 



(47) 



It is convenient to introduce at this point the symbol J\f(s,t) defined in Eq. @. It allows 
to rewrite (f|7|) with the help of some standard transformations of binomial coefficients as 



10 



n-1 



K 2 (n;n/4) 



~^ on+l 



9=1 

n-1 



n 



q 



M(q, n — q — 1) 



+ E Mfc n - g) + (-l)"JV(n - g, g)]' 

9=1 



(48) 



This expression is displayed in Fig. [I] together with (|4T) ) in order to illustrate the equivalence 
of the two results for the form factor. 



1.5 




1 1 1 1 1 

10 20 30 40 

n 

FIG. 1. Form factor for the ring graph (see inset). The crosses and the connecting heavy full 
line show the two equivalent exact results (|40| ) and (|47| ) for rj = it/ 4. The thin dashed lines represent 
the approximation ([Tl|) . The heavy dashed line exhibits the form factor of a CUE ensemble of 2 x 2 
random matrices (f43|), which can be obtained by an appropriate averaging over rj (see text). 

An independent proof for the equivalence of ([TOP, can be given by comparing the 
generating functions of K 2 (n]7T/4) in the two representations ||. We define 

oo 

G(x) = K 2(n; tt/4) (2x) n (\x\ < 1/2) (49) 

x=l 

and find from fl4*0|) 

l-2x 2 ^ 2 \mj 
1 l-2x 11 -6x 
" 2 VI + 4x 2 ~ 2 1 - 2x ' ( ^ 
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On the other hand we have from (|48| 



G(x) 



x 



1 — X 



G 1 (x) + G 2 (-x) 



(51) 



with Gi^{x) defined in the introduction. A convenient starting point to obtain the r.h.s. of 
(0), (HI) is the integral representation 



Af(s,t) 



-lY 



2vri 



dz(i + z-ya- z 



\s-l 



(52) 



where the contour encircles the origin. With the help of (|52"D we find 

oo 

g(x,y) = £ N(s,t)x s y t 

s,t=l 

■1 oo « oo 



s,t=l 



:,t=l 



2vri 



s,t=0 



7 V idZ" 77" 

i„4±J 1-xl 



1 + 2 



(1-z) z + y(l + z) 



xy 



+ -x + y-2xy) 



(\x\,\y\<i/V2) 



(53) 



which was already stated in @. The contour \l + z 1 \ = |1 



v^2 has been chosen such 



that both geometric series converge everywhere on it. Now we have 



1 r c\ 7 f\ 2-' 

Gl(l2) = WJ~^ ^ N{ S ,t)N{s',t'){xzy{x/zr {xz'f^/zj 



(2vri) 

X 4 

(2vri) 5 



,t=l s',t'=l 



dz dz' 



(1 + X2')(l + x[z' - z] - 2x 2 zz') (z' + x)(zz' + x[z - z'\ - 2x 2 ) 



, (54) 



where \x\ < 1/ \/2 and the contour for z, z' is the unit circle. We perform the double integral 
using the residua inside the contour and obtain @ and in complete analogy also (f|) such 
that finally 



G{x) 



x 



+ 



X 



1 Ax 2 - 2x + 1 



l-x 2x-l \^4x 2 + l l-xj 2(l-2x) v / 4x 2 + l 2 



(55) 



The proof is completed by a straight forward verification of the equivalence between the 
rational functions (^) and (|55]). 

The identities (||), (0) follow now by separating even and odd powers of n in ([5]) and 
(fl7|). In terms of Kravtchouk polynomials these identities can be written as 



2m- 1 

E 



/2m - 


;) 


\2m - 





2m 



-P 



2m— 1,2m— q 



(?) 



i2m+l 



+ (-!)' 



2m 



m 



- 2 



(56) 



and 
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^/ 2m ^ 



2m + 1 

q 



P, 



2m,2m+l—q 



(?) 



>2m+2 



-2(-l) 



m / 



(57) 



Finally we will derive the CUE result ( fi3"D for the ensemble of graphs defined in the previous 
subsection starting from the periodic-orbit expansion (|4]). We find 



rn/2 

K 2 (n) = dfJ,(r))K 2 {n;r)) . 
Jo 

Inserting (f44|), expanding into a double sum and using 



(58) 



tt/2 



drj sin 



2{v+v')+l 



7] COS 



2(n-i/-i/)+l 



V 
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2(n + 1) \v + v' 



(59) 



we get 
# 2 (n) 



n + 1 



(60) 



n-l 



4(" + 1) £ 5 



(-1)^' 



z^z/ 



'g — l\/n — g — l\/g — l\/n — g — 1 
.iz-ljl JU'-lJl z/-l 



Comparing this to the equivalent result we were led to the identity ([!]) involving a 
multiple sum over binomial coefficients. In this case, an independent computer-generated 
proof was found [Til, which is based on the recursion relation 



q 2 F uy (n, q )-(n-q- l) 2 F v y(n, q + 1) + (n - l)(n - 2g - + 1, g + 1) = . (61) 



This recursion relation was obtained with the help of a Mathematica routine , but it can 
be checked manually in a straight forward calculation. By summing (|6ll) over the indices 
v, u', the same recursion relation is shown to be valid for S(n, q) |TT|Jr2]] and the proof is 
completed by demonstrating the validity of (|l|) for a few initial values. Having proven ([[]) 



we can use it to perform the summation over u, v' in (|6(]) and find 

-i n— 1 



K 2 (n) 



n 



1 9 =1 



n 



+ 



n 



rr 



n + 1 n + 1 



(62) 



which is now obviously equivalent to the random-matrix form factor 



D. Trace identities 

Let S be an arbitrary unitary matrix with a non degenerate spectrum and s n = trS r 
Then 



lim e V s* s n+u e n€ = s u , 

£—♦0 

n=no 



(63) 



for arbitrary integers Uq and [13 
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We shall now apply this identity to the ring graph with rj = tt/4 in order to prove @. 
Again, the traces of S n can be expanded in periodic orbits which are grouped into n + 1 
families with equal phases $(n, q) = q<pi + (n — g)0 2 (0 < <7 < w). Thus, one can write 

s n {k) = j2MnA)^ q \ (64) 

9=0 

where A(n, q) is the coherent sum of all amplitudes of PO's in the corresponding set. A(n, q) 
can be expressed in terms of Kravtchouk Polynomials as 

1 | 1 for q = or n 

A(n,g) = 7F\(-irn^)(r;) 1/2 ^-^) fo r o< g <n (65) 



(compare Eq. (|5|)). This is equivalent to (|8|). Substituting (|64j) into the trace identity (|63D, 
we get for arbitrary integers v and Uq 



oo n n+v 



lim e £ e~ ne J2J2 A (n + v,p)A(n,q)e 



i[{p-q)4>i+(v-{p-q))<l>2\ 



u g=0 p=0 



= A & K ) e i[ ^ 1+(i/ - K)</,2] . (66) 

K=0 

This is valid for arbitrary phases 0i and 02 and therefore the coefficients of the phase factors 
e i*(f>K) on both sides are equal. Eq. (|9|) follows. 



IV. CONCLUSIONS 

We have shown how within periodic-orbit theory the problem of finding the form factor 
(the spectral two-point correlation function) for a quantum graph can be reduced exactly 
to a well-defined combinatorial problem. Even for the very simple graph model that we 
considered in the last section the combinatorial problems involved were highly non-trivial. 
In fact we encountered previously unknown identities which we could not have obtained 
if it were not for the independent method of computing the form factor directly from the 
spectrum. However, since the pioneering work documented in [|12| the investigation of sums 
of the type we encountered in this paper is a rapidly developing subject, and it can be 
expected that finding identities like (0), (|J) and (0) will shortly be a matter of computer 
power. 

Numerical simulations in which the form factor was computed for fully connected graphs 
(vi = V — 1 V?) and, e. g., Neumann boundary conditions indicate that the spectral corre- 
lations of 5* match very well Dyson's predictions for the circular ensembles COE or CUE, 
respectively 0. The agreement is improved as V increases. We conjecture that this might 
be rigorously substantiated by asymptotic combinatorial theory. A first step towards this 
goal was taken in the present paper. 
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